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the hydrodynamic and thermodynamic character within the ejector. The paper first briefly describes ejector 
including fundamental principle, flowing and mixing mechanism and the method of model establishment. 
Then various models consisting of ideal assumptions, governing equations, auxiliary conditions, solution 
methods and main results are presented. The models can be classified into two main categories: (i) steady 
thermodynamic models which can be further subdivided into single-phase flow model and two-phase flow 
model and (ii) dynamic models which are also subdivided according to the flowing phases considered. It has 
been shown that the dynamic models have higher prediction precision and give more information 
compared with the steady thermodynamic models. In addition, the simplified empirical and semi-empirical 
models based on measured data are briefly discussed. This review is useful for understanding the evolution 
process and the current status of the mathematical models on ejector and highlighting the key aspects of 
model improvement such as the mixing mechanism, the capture of the shock wave, etc. 
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Nomenclature 


cross-section area (m7) 

local sonic velocity (m/s) 
constant-pressure heat capacity (kJ/kg K) 
diameter (m) 

total energy, Eq. (121) 

friction factor, Eqs. (42), (45), (46), (92) 
specific enthalpy (kJ/kg) 

turbulent kinetic energy, Eqs. (4), (124) 
length (m), Eqs. (45), (46), (92) 

Mach number 

mass flow rate (kg/s) 

pressure (MPa) 

condensation rate, Eqs. (108), (112) 
Reynolds number, Eq. (42) 

universal gas constant (kJ/kg K) 

entropy (kJ/kg) 

temperature (K) 

velocity (m/s) 

velocities in various directions (m/s) 
location along the ejector axis(m), Eqs. (4), 
(94-96); quality of the fluid (kg/kg) 
axial coordinate in cylindrical geometry (m) 


Greek letters 


é 


WE HC ADTES YON 


void fraction, Eqs. (106), (109), (110), (115), (117); 
rate of turbulent energy dissipation 
heat capacity ratio 

dissipation energy (kJ/m? s), Eq. (136) 
production rate 

isentropic efficiency 

dynamic viscosity (Pa s) 

density (kg/m?) 

shear stress 

specific volume (m?/kg) 

area contraction ratio 

mass flow ratio 

momentum loss factor, Eq. (92) 


Subscripts 


c 


d 
e 
f 
g 


condenser 

diffuser 

exit; evaporator 

saturated liquid, Eqs. (32) and (33) 
generator; saturated vapor, Eqs. (32), (33); 
gas phase 

inlet, location along ejector axis 
isentropic process, Eqs. (43), (44) 
liquid phase 

phase change 

mixture; mixing section 

primary nozzle 

nozzle throat 

primary stream 

secondary stream 

saturation 

conditions downstream of shock wave 


sub sub-cooling 

t total temperature/pressure 
v vapor 

y mixing section 


1, 2,3,4 locations along the ejector 


1. Introduction 


The ejector refrigeration system has great market potential 
because it offers some remarkable advantages which include: (i) it 
can alleviate environment problems by using low grade thermal 
energy sources such as solar energy, geothermal energy and waste 
heat to drive the system instead of high grade electric energy, 
hence it can reduce CO; emissions resulting from the combustion 
of fossil fuels, (ii) it is simple and with no-moving parts, which 
makes it noise-free, reliable, long lifetime, low initial and running 
cost as well as require practically no maintenance and (iii) natural 
substances such as water can be utilized as working fluids, which 
have zero ozone depletion potential. The main flaws of the system 
are low performance efficiency and rigid requirements on the 
operation parameters which limits their widespread application. 
The most direct way to solve this problem is to investigate the 
ejector, which is the most critical component where highly 
irreversible transformations take place. In particular, due to its 
importance on the whole system, the ejector can be taken as the 
core of an optimization analysis aiming to define its operational 
conditions and the geometric configuration that make the COP of 
the system maximum [1]. Some review papers have been 
published to summarize the research efforts and achievements 
focused on ejector performance optimization. Chunnanond and 
Aphornratana [2] summarized the ejectors and their applications 
in refrigeration. They concluded that the understanding of the 
ejector theory had not been completely cleared. In order to 
enhance the efficiencies and reduce the cost of ejector cooling 
systems, efforts made by several researchers have been summar- 
ized in yet another review paper [3]. The review by Sun and Eames 
[4] outlined the developments in mathematical modeling and 
design of jet ejectors. However, the review considered only the 
thermodynamic model based on two basic approaches, the mixing 
of the primary fluid and the entrained fluid either at constant- 
pressure or at constant-area. 

With the development of computer hardware and numerical 
methodology, mathematical models are being used for better 
understanding of the compression process, system design and 
performance evaluation considering the hydrodynamic perfor- 
mance, condensation and two-phase flow. The advantages of this 
method are that it takes less time and cost than experimental 
method for predicting the performance of an ejector. The second 
point is that mathematical models can produce large volumes of 
results at virtually no added expense and it is very convenient to 
perform parametric research and optimization analysis. The third 
reason is that by experimental means, some parameters are 
difficult to be obtained. The last but not the least reason is the 
fundamental physics of ejector used in various fields are similar, so 
that other installations may be investigated too by utilizing the 
model. Consequently, constructing valid mathematical models of 
the ejector has become the key subject of many studies. Lots of 
mathematical models have been constructed and employed to 
analyze, develop and design ejectors. This article examines the 
progress made in the area of mathematical modeling on the 
ejector, and summarizes comprehensively the numerous signifi- 
cant works that has been done on modeling the ejector. 
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Based on the conservation principles of energy, mass and 
momentum, the hydrodynamic and thermodynamic performance 
of ejector can be formulated with complex mathematical detail as 
governing equations and auxiliary relations. This is the main and 
universal method for constructing mathematical model. There are 
also some simplified models for ease of computation such as semi- 
empirical model [5]. Moreover, it is an effective way to get the 
empirical performance fitted correlation using large amounts of 
experimental results. But its application is restricted in the specific 
range of experimental working conditions. This paper focuses on 
the review of the mainstream mathematical models. Besides, some 
empirical/semi-empirical models are briefly introduced. 

The first section of the paper details model establishment basis, 
which can be either thermodynamics or dynamics arranged 
according to the flowing phases. For each model reviewed, attempt 
has been made to gather the assumptions, governing equations, 
auxiliary relations and the solution techniques for ensuring the 
integrality of the model. The results and validation of some models 
are listed where comparison between them are needed. 


2. Problem description 
2.1. The fundamental principle 


An ejector is also known as jet, injector or jet pump in different 
literatures. The main components of an ejector include a primary 
nozzle (also named as motive nozzle in some literatures), the 
suction chamber, the mixing chamber (including a convergent 
chamber if available and a constant-area throat tube which is 
named as secondary throat in some literatures) and the diffuser 
(Fig. 1). 

The primary nozzle usually is a convergent-divergent nozzle. 
As the high pressure fluid, known as “a primary fluid” or “motive 
fluid”, expands and accelerates through the primary nozzle, it 
flows out with supersonic speed to create a very low pressure 
region at the nozzle exit plane and subsequently in the mixing 
chamber. Hence, a pressure difference between the streams at the 
nozzle exit plane and the secondary fluid inlet is established and “a 
secondary fluid” or “entrained fluid” is drawn by the entrainment 
effect. 

Based on Munday and Bagster’s theory [6], it is assumed that 
the primary fluid flows out without mixing with the secondary 
fluid immediately. It expands and forms a converging duct for the 
secondary fluid. At some cross-section along this duct, the speed of 
the secondary fluid increases to sonic value and chokes, were 
named as “hypothetical throat” (or “fictive throat” or “aero- 
dynamic throat” in some literatures). Then the mixing process 
begins after the secondary flow chokes (double-choking, critical- 
mode operation). It is assumed that the pressure of the two 
streams is uniform at the mixing section. Also, this mixing causes 


Nozzle Mixing 
cti. section 
a section rls 


Constant-area 
section >a 


the velocity of the primary flow to be stepped down while the 
secondary flow is accelerated. By the end of the mixing chamber, 
the two streams are completely mixed and the static pressure is 
assumed to remain constant until it reaches the constant-area tube 
section. Due to a high pressure region downstream of the mixing 
chamber’s throat, the flow undergoes a succession of normal and/ 
or oblique shock waves, also called a shock train, accompanied by a 
corresponding pressure rise. This shock causes a major compres- 
sion effect and a sudden drop in the flow speed from supersonic to 
subsonic. A further compression of the flow is achieved as it is 
brought to stagnation through a subsonic diffuser. 

It can be observed that due to the special geometrical structure 
of the ejector, the flow in the ejector is very complex, which 
involves supersonic flow, shock interactions, and turbulent mixing 
of two streams and two-phase flow in some cases. Hence the main 
endeavor of modeling the ejector is to reasonably explain the 
complex flow and mixing mechanism. 


2.2. Flow and mixing mechanism in ejector 


As mentioned in Section 2.1, the flow in the ejector can be 
divided into three sections: 


(1) Nozzle section: A primary fluid is accelerated to supersonic 
speed by the convergent-divergent primary nozzle, which 
forms low pressure region at the nozzle exit plane and 
produces the entrainment effect to entrain the secondary fluid. 
Since the velocity is very high, the process can be considered as 
adiabatic and the isentropic relations are applicable. 

(2) Mixing section: Two feasible mathematical methods, the 
constant-area mixing model and constant-pressure mixing 
model proposed by Keenan et al. [7] are used to describe 
the mixing mechanism. Also, a normal shock wave may 
occur if a supersonic fluid mixture enters the diffuser. In this 
case, a sudden reaction in the mixture velocity and a rise 
in the pressure take place. In addition, the fluid mixture 
easily undergoes phase change and a condensation shock 
may occur. 

(3) Diffuser section: The mixture of primary and secondary 
flows passes through the diffuser, and converts kinetic 
energy into pressure energy. At the diffuser exit, the velocity 
is reduced to zero and the pressure is lifted high enough to 
cause discharge. 


Therefore, the mathematical description of the flow inside the 
ejector is very complex. Besides the conservation equations of 
mass, energy and momentum, some gas dynamic equations, state 
equations, isentropic relations as well as some appropriate 
assumptions need to be used to assist in the description of the 
flow and mixing in the ejector. 


Diffuser section 


Primary fluid a 


Secondary fluid 


Shock 


Mixed flow 


Fig. 1. Principle structure of ejector. 
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2.3. Mathematical model establishment 


The main steps used to establish mathematical models include: 
(i) assumption proposal, (ii) deriving the governing equations 
based on the mass, energy and momentum balances, (iii) providing 
auxiliary relation to close the governing equations, and (iv) 
adopting proper mathematic approach to solve the model. For 
clarity, the system is considered to be one-dimensional (for 
thermodynamic models). Explanations of the main terms in the 
governing equations are presented below. 

Some ideal assumptions are needed because the choking, shock 
and mixing of the two streams occurring in the ejector are too 
complicated to understand completely. On the other hand, some 
factors that do not influence the flow significantly can be 
neglected. Then the complexity of the governing equations can 
be reduced and the solution is relatively simple and time-saving. 
Some basic assumptions made before establishing mathematical 
models are as follows: 


(1) The inner wall of the ejector is adiabatic. 

(2) The flow inside the ejector is steady and isentropic. 

(3) The primary fluid and secondary fluid are supplied to the 
ejector at zero velocity. 

(4) The velocity at the ejector outlet is neglected. 

(5) The two fluids begin to mix with a uniform pressure at the 
mixing section. 


From assumption (1), the heat transfer between the ejector and 
the environment can be neglected and the isentropic expressions 
can be used based on assumption (2). Whereas from assumption 
(3), the inlet fluid is at a stagnant state, therefore, the inlet pressure 
and temperature of the fluid are also equal to the total pressure and 
temperature as well. From assumption (4), the fluid is brought 
down to stagnation when it discharges out of the ejector. 


2.3.1. General governing equations 
Under ideal condition, the general governing equations for the 
nozzle, mixing and diffuser sections are as follows: 


5 PUA; = 5 PelleAe (1) 


Conservation of mass : 


Conservation of momentum: PA; + 5 miui 
= PeAe + X Meule (2) 
Conservation of energy : 5 mi(hi + u? /2) 
= $ me(he + ue/2) (3) 


2.3.2. Auxiliary conditions and mathematical methods 

In order to close the equations, some auxiliary conditions have 
to be provided. They are boundary conditions, initial conditions for 
iteration, gas dynamic equations, entropy equation and the state 
equation. Also, the thermodynamic and transportation properties 
of the fluid should be given as well. 


(1) Boundary conditions: Boundary conditions describe the 
behavior of the simulation at the edges of the simulation 
region. The boundary conditions commonly encountered in 
ejector models, are usually in pressures at inlet and exit of the 
ejector. Also the mass flow rates or the velocities of the primary 
and secondary fluids have been used as boundary conditions in 
some literatures. 

(2) Initial conditions: The initial condition used in ejector models 
is the expansion ratio. The entrainment ratio and/or the cross- 


section area of the constant-area throat tube are also used for 
global models. In addition, the inlet conditions at the boundary 
mesh were used for local models. 

(3) Turbulence modeling: As the flow in the ejector is from 
subsonic to sonic and then to supersonic, the turbulence has 
strong effect on the flowing process. Generally, for dynamic 
models, the Boussinesq hypothesis is used for turbulence 
(Eq. (4)). This means that they are based on an eddy viscosity 
assumption, which makes the Reynolds stress tensor resulting 
from equation averaging to be proportional to the mean 
deformation rate tensor. 


Ou; _ Ou; 2 l Ou; 
Hila | ax D 3 (pk F Ht Sx Ou (4) 


puu, 


For thermodynamic models, the turbulence cannot be 
modeled in detail. The dissipation term is interpreted by 
frictional and mixing losses which are taken into account by 
using corresponding coefficients introduced in the governing 
equations. Generally, these coefficients need to be determined 
experimentally. 

(4) Auxiliary relations: 


Definition of Mach number: M = u/c (5) 


Definition of sonic velocity: c 


=4/ 6) , for ideal gas, c 
= YR (6) 


Based on isentropic assumption, the following gas dynamic 
equations are frequently used in the thermodynamic models: 


(vi+1)/(yi—1) 
m PVT y 1 ae (7) 
Ai JReTe N UFM? 
y+1 
Air Mi [1+4 Miu |0 (8) 
A Mia | 1+%51M’? 
P 1+5 Mia \ 
i E eae (9) 
Piyt 14+55M; 
Tr -Ylin 
Ta 1+ ~y Mu (10) 
: P 
Entropy equation : o = constant (11) 
State equation: h = h(P, p) (12) 
For constant pressure mixing, h= CpT (13) 


(5) Mathematical solution: Due to the complexity and difficulty in 
modeling the flow in the ejector completely, most of the 
mathematical models are one-dimensional thermodynamic 
models which only calculate the steady-state explicit equa- 
tions and obtains state and operation parameters along the 
ejector. Detailed information such as shock interactions, 
turbulent mixing of two streams cannot be obtained by these 
models. With the development of mathematical method and 
computer science, especially with the birth of commercial CFD 
software, numerical methods have been applied to solve the 
partial implicit differential governing equations which can give 
more information. At present, the most widely used solution 
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method for ejector model is finite difference method. The 
governing equations are solved using the commercial CFD 
software package FLUENT with the governing equations being 
discretized using a control volume technique. 


3. Thermodynamic model 


Among the different models, the model based on thermo- 
dynamic analysis is the simplest, which solve the equations in one 
dimension. In this model, the isentropic relations and some gas 
dynamic equations are used to assist in the description of the 
relationship between the temperature, pressure, enthalpy and 
Mach number and/or velocity. Hence, most of them are explicit 
steady-state models and work at design conditions. 

Moreover, in order to simplify modeling, most of the models 
assume that the inlet fluid is superheated and no phase change 
occurring in the ejector. These models are developed based on 
single-phase flow. However, the experimental results indicate that 
the fluid in the ejector easily undergoes condensation shock. Also, 
for some practical application, the ejector works as a pump by 
using high pressure water to entrain gas or high pressure vapor to 
entrain water. In these cases, the flow in the ejector is two-phase 
flow. Therefore, in the next section, the models have been 
subdivided into single-phase flow models and two-phase flow 
models. 


3.1. Single-phase flow model 


3.1.1. Constant-pressure mixing model 

Since the birth of the first steam jet refrigeration system, its 
advantages have made it to enjoy popularity. In 1942, the 
theoretical prediction of the ejector progressed further, when 
Keenan and Neumann [8] established a one-dimensional con- 
tinuity, momentum and energy equations to predict the perfor- 
mance of the ejector. However, the difficulty in offering a 
reasonable analytical solution for the momentum conservation 
equation (Eq. (14)) during mixing, made the use some experi- 
mental coefficient necessary in model. 


DF = $ AdP = m9(1 | @) Um 


Therefore, in Keenan et al.’s following work [7], they introduced 
two feasible theoretical methods to solve the problem, the 
constant-pressure mixing model and the constant-area mixing 
model. The former considers that dP = 0 during mixing while the 
latter considers dA = 0. 

Keenan et al.’s work lays a foundation for one-dimensional 
design theory of the ejector. They pointed out that a constant- 
pressure mixing ejector gives better performance than a constant- 
area mixing ejector. Therefore, most of the mathematical models 
that followed are based on the constant-pressure mixing ejector. 
However, Keenan et al.’s initial model could not give detailed 
information on the choking phenomena that is easy occurring in 
the supersonic ejector when the back pressure is low. 

Munday and Bagster [6] further developed the constant- 
pressure mixing model by assuming that the primary fluid flows 
out without mixing with the secondary fluid immediately after 
discharging from the nozzle exit. They assumed that both fluids 
mix somewhere at the downstream of the nozzle exit (Fig. 1, 
section y-y) in the suction chamber and a “hypothetical throat” for 
entrained fluid is formed here. 

However, neither Keenan’s nor Munday’s model has taken 
consideration of irreversibility due to friction. Eames et al. [9] 
modified Keenan’s model to include irreversibility associated 
with the primary nozzle, mixing chamber and diffuser based 


MpU2p — MsUgs (14) 


on constant-pressure mixing theory but without considering 
the choking of the secondary flow. For simplification, besides 
assumptions (1-5) mentioned in Section 2.3, they further 
assumed: 


(1) Friction losses were introduced by applying isentropic 
efficiencies to the primary nozzle, diffuser and mixing 
chamber. 

(2) The two streams mixed at the primary nozzle exit plane. 

(3) Mixing of the two streams is complete before a normal shock 
wave can occur at the end of the mixing chamber. 

(4) The working fluid is ideal gas. 

(5) The flow is one-dimensional. 

(6) The constant-pressure mixing occurred at the inlet of the 
constant-area section (Fig. 2). 


Based on these assumptions, the model was then written as: 
Primary nozzle: 


2 


u 
Energy equation : = = Na (hpt — Mp1) (15) 


Mp1 = oo BY | (16) 
mora) a7) 


where hy, is the flow enthalpy at the nozzle exit for isentropic 
expansion. 

Mixing section: Based on assumption (6), P4 = Pm, A1 = Am, and 
assumption (5), the momentum equation (Eq. (2)) between section 
1-1 and section m-m can be simplified into Eq. (18) with nm as the 
efficiency for the whole mixing chamber. 


Nm (MpUp1 + MsUs1) = (Mp + Ms)Um (18) 


Using Eq. (6), Eq. (18) then can be written in terms of the Mach 
number: 


M: — Mp1 + OMS; VTs/Tp 
m (@+ 1)(14 @T;/Tp) 


(19) 


where œ is the entrainment ratio: w=m,/mp, and 


M* = ye M?)/(1 , ut M’), it is the actual mixture velocity 


divided by the sonic velocity of the mixture at critical conditions. 

Shock wave: The relations between the Mach number 
and pressure lift upstream and downstream of the shock 
wave are 


2 2 
= _ 2 Ff _ M2 _ 
Ms = (25 +ma)/ (PMs 1) (20) 
14 yM2 
Ps/Pm =TM (21) 


Diffuser section: The pressure lift in the diffuser described in 
Mach number including nq as efficiency is 


aM, 
P4/P3 = (m È 1)’ (22) 
Give the temperature, pressure and mass flow-rate of the 


primary and secondary fluids and assuming a value of P1/Pst, 
calculating Eqs. (16), (17), (19-22) consequently, then you can 
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Suction 


Nozzle 


section 


Primary fluid 


Secondary fluid 


Constant-area 
section 


Diffuser section 


Mixed flow 


Fig. 2. Ejector model used by Eames et al. 


obtain the exhaust pressure P4, repeat the calculations with a new 
value of P,/P,+ until the maximum P; is obtained. The primary 
nozzle, diffuser, and mixing chamber efficiencies of 0.85, 0.85 and 
0.95, respectively were used for the analysis. The authors verified 
the theoretical results with experimental data and found that 
choking of the secondary flow in the mixing chamber of the ejector 
plays an important role in the system performance. Maximum COP 
can be obtained when the ejector is operated at its critical flow 
condition. 

Based on the theory proposed by Munday and Bagster [6] 
and using the model developed by Eames et al. [9], Aly et al. 
[10] presented two models. One was developed by applying 
steady-state equations of energy, momentum and continuity at 
the nozzle, diffuser and mixing section to determine the 
pressure and velocity at each section from which the system 
performance can be obtained. Based on the assumptions (1-5) 
mentioned in Section 2.3 and taking into account of the losses 
due to wall friction in the nozzle, and those in the mixing and 
diffuser sections, the mathematical description of the model is 
as follows: 

Primary nozzle: 


= Nn (hp, — Mp1) (15) 


Mixing section: In the mixing section, both the momentum and 
energy balance equations were used. Comparing with Eq. (18), the 
momentum term of the secondary fluid at left side is omitted 
(Eq. (23)) based on assumption (3) in Section 2.3. The energy 
balance is shown in Eq. (24). 


NmMpU = (Mp + Ms)u2 (23) 


u3 uy 
(Mp 4 ms) (a j 2) = msh; 4 p(y H 3) (24) 


Shock wave: The Eq. (25) similar to Eq. (20) was derived to 
calculate the mixture velocity after the shock wave 


tig (Fae) Pam) xc (25) 


By applying the energy and continuity balance before 
and after shock wave, the enthalpy and density after the shock 
wave were 


u? už 
h; =h +4 -3 (26) 


u 
P3 = p= (27) 
u3 


Diffuser section: The energy equation taking account the 
isentropic efficiency for the diffuser section is 


TE E (28) 
AS 2na 
where hai; is the flow enthalpy for the isentropic process. 

Assuming a certain value for the entrainment ratio and 
considering the isentropic expansion process in the primary 
nozzle, the pressure P; and enthalpy h, at the nozzle exit can be 
obtained. Then Eqs. (15), (23-28) can be used to obtain the diffuser 
exit enthalpy h4. The calculation is repeated with a new 
entrainment ratio until the required exit pressure is obtained. 

The second model used the assumptions in the first model and 
in addition, the flow inside the ejector was considered as ideal gas 
with constant specific heat ratio both for the superheated region 
and the wet region separately. Following Eames et al.’s model of a 
steam-vapor ejector [9], they calculate the pressure ratio and 
Mach number of the flow in the ejector to study the characteristics 
of the system. 

For validity, the results of the two models were compared with 
the empirical correlation of Power [11] and good agreement was 
found. Then both models were used to analysis the effect of the 
efficiencies such as nozzle and diffuser efficiencies, as well as the 
design parameters on the performance of the ejector. They found 
that the nozzle and diffuser efficiencies have a significant influence 
on the system performance. 

In order to take account of the chocking of secondary fluid that 
neither Eames et al’s nor Aly et al’s models considered, Huang et al. 
[12] presented a critical-mode (double-chocking) model based on 
Munday and Bagster’s theory [6] by assuming that the primary 
fluid mixed at some cross section (Fig. 2, section y-y) downstream 
of the nozzle exit. Hence, the model is more complex, besides all 
equations derived by Eames et al., the gas dynamic relations 
(Eqs. (6-10)) and state equation Eq. (13) were used in the 
calculation. In addition, the area relationship (Eq. (29)) and the 
energy conservation in mixing section (Eq. (30)) were used as well. 


Apy + Asy = A3 (29) 
u2 u2 
Mp (cot + P) +m; (cots + 7) 
2 
= (mp + ms) (CoTm cf a) (30) 


Eleven different ejectors were tested to determine the 
coefficients used in the model and to verify the theoretical 
analysis. Then the performance analysis was carried out to 
determine the entrainment ratio and the required cross sectional 
area of the constant-area throat tube. In this model, the initial 
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Fig. 3. Ejector and the Molier chart of an ammonia ejector. 


conditions were the throat and exit area of the primary nozzle, the 
stagnant temperature and pressure of the primary and secondary 
fluid and the critical back pressure. 

It should be pointed out that the previously mentioned models 
are based on ideal gas assumption which does not reflect the actual 
process occurring in the ejector. To solve this problem, Rogdakis 
and Alexis [13] improved the model proposed by Munday and 
Bagster [6] by using the thermodynamic and transportation 
properties of real gas. The model was based on the Molier chart 
of an ammonia ejector (Fig. 3(b)). The momentum conservation 
during mixing was employed in the model. Detailed governing 
equations are as follows: 

Primary nozzle: 


The energy balance: u2,/2 = hy — hat (31) 


The state equation: $s, = s(T1, P1) = Sas 
= XasSgas + (1 — Xas) Stas (32) 
and 


h= Xashgas + (1 — Xas)hfas (33) 


The isentropic efficiency is 44, = (hı — ha1)/(hı — has) (34) 


Similar equations were developed for the secondary fluid. 
Mixing section: 


Momentum balance: wai + @Ua2 = (1 + œ)up (35) 


Similarly, by employing the energy balance and state equations 
Eqs. (31) and (32), the enthalpy h,,hg and the specific volume ve, Va 
can be obtained. 


Diffuser section: 


sa = S(Ta, Pa) (36) 
Sas = Sa = S(T3s, Pa) (37) 
həs = h(T3s, Pa) (38) 
nas = (ha — h3s)/(ha — hs) (39) 
h3 = h(T3, P3) (40) 


Given the pressure for the state (1-3) (Fig. 3(b)) and using 
iterative calculation, the maximum entrainment ratio wœ can be 
found with constant superheated temperature T4. 

The ejector efficiency of compression is: 


N = (Mg + Me) (h3 — h2)/(Mg(hı — has)) 
= (1 + œ)(h3 — h2)/(hı — has) (41) 


The subscripts 1, 2, 3, b, c, d, a1, a2, as and 3s are the cycle 
locations as illustrated in Fig. 3. 

The model can be used for one component or two components 
ejector. They used the model to analyse the performance of an 
ejector cooling system providing air conditioning and it was found 
that the shock phenomena on the entrained vapor play important 
role in ejector performance and the COP of the system. Also, from 
their work, they observed that the highest operating efficiency can 
be achieved if the generator pressure increases with decreasing 
condenser temperature, when the evaporator temperature is fixed. 
In addition, they showed that the maximum coefficient of 
performance is a linear function of generator temperature, a 
quadratic function of condenser temperature and a cubic function 
of evaporator temperature. 
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Fig. 4. Ejector model used by Selvaraju and Mani. 


The friction loss is closely related with the velocity, and the 
velocity varies in a big range along the ejector. Therefore, using a 
constant coefficient to represent the friction loss as in above- 
mentioned models is a simplified treatment. Also, Aly et al. [10] 
found that the efficiencies have a significant influence on the 
system performance. Taking this into account, Selvaraju and Mani 
[14] developed a model based on Munday and Bagster’s [6] theory 
for critical performance analysis of the ejector system (Fig. 4) by 
applying an expression to describe the friction loss in the constant- 
area section: 


1 
= 2.0log(Rem,/ fm) — 0.8 (42) 
Vim oa 
Then they derived the minimum area per mass flow rate and the 
velocity for both primary and secondary fluid at the primary nozzle 
outlet and the aerodynamic throat, respectively. 


A 


<=? (43) 
Mr /2np(Ihg — hi)is 
Upe = 4/2Np(ħg — Ape). (44) 
The velocity of the mixed fluid is expressed as, 
Um = (MpUpe + MsUme) + (Pse — Pm)Am (45) 


mm (1 +p) 


The entrainment ratio and the characteristic area ratio for 
ejector are, 


Upe — War — 3Uar fm(lm/Dm) 


w= 46 
Uqi — Use + dud fm(lm/Dm) ! i 

Am UmUt 

— = (1+ 0) —— 47 

a ( as (47) 


The model was used to analyse the effect of the compression 
ratio and the driving pressure ratio on the critical entrainment 
ratio and to predict the performance of an ejector refrigeration 
system using various refrigerants. They found that the entrainment 
ratio is reversely proportional to the compression ratio, propor- 
tional with the driving pressure ratio, and the system with R134a 
gives better performance. 

By using a model similar to that of Rogdakis and Alexis [13], Yu 
et al. [15] derived a thermodynamic model based on isentropic 
process to estimate the entrainment ratio of the ejector for both a 
refrigerator and a pump: 


oo = Yamal (Apt — Mpt)/(Hmas — mm) ~ 1 (48) 


where nnm and na are the efficiencies of the nozzle, mixing 
chamber and the diffuser. They can be obtained from the following 
expressions: 


Apt = hp1 
_ 49 
In hpt = hp1 is ) 
Um 
Nm = up 50) 
hma s= hmm 
_ "ma, 51 
a hma = hmm ) 


This model was used to predict the ejector refrigeration system 
with mechanical sub cooling which uses an auxiliary liquid-gas 
ejector to enhance sub cooling for the refrigerant from condenser. 
It was found that the performance of this cycle improved more 
than the traditional ejector refrigeration system and an optimum 
mechanical sub-cooling degree exists for different refrigerant and 
operating conditions. 

In order to obtain more detailed information on the ejector, 
Ouzzane and Aidoun [16] derived a local mathematical model 
and computer programs for ejector studies in refrigeration 
cycles, one program for optimal ejector design and the other for 
simulation with more in-built flexibility. The model is based on 
Munday and Bagster’s theory [6] and isentropic flow in the 
nozzles and the diffuser. The control volume along the ejector 
axis is shown in Fig. 5. The conservation equations for the 
control volume are, 


Energy conservation: h; + i =h_, +5u (52) 
Mass conservation : p,UjAj = P;_1Ui_1Ai-1 (53) 
Isentropic condition: Si = $;_1 (54) 
Momentum conservation: ,u? + P; = pj 1U? 1 + Pii (55) 


The iteration variation for the solution procedure is Mach 
number, such that 


Mi = M,_, + AM (56) 
\ 
| 
Mit | | M; 
hi | | n 
Pii l P; 
w | ly, 
a E 
Pii | | pi 
a E 
aia | T 
Aia l | A 


Fig. 5. Control volume used in Ouzzane and Aidoun’s model. 


1768 S. He et al./Renewable and Sustainable Energy Reviews 13 (2009) 1760-1780 


The enthalpy value obtained from Eq. (52) is corrected for 
friction as: 
hi = hy_y + (his — His) (57) 

Together with the defining of Mach number and sonic velocity, 
and given the initial condition, the model can be calculated step by 
step. 

The models were developed for qualitative and quantitative 
studies on the operation and performance for optimal and off 
design operation. The model can also be used to represent the 
global characteristics and detailed evolutions of all the relevant 
parameters. In addition, it can be used to study the influence of 
individual parameters on the general ejector operation. They used 
the model to analyze the effects of geometric configurations on 
ejector operation under different conditions showed that increas- 
ing ejector diameter improved entrainment ratios and capacities 
but exit pressures decreased accordingly. Component length 
affects not only the overall ejector size but also the capacity and 
the stream mixing process. The mixing chamber length was 
particularly important for the control of the shock wave intensity. 
Less irreversibility and limited superheat was produced in the 
ejector when the length was adjusted to bring supersonic mixed 
flow to near sonic conditions for maximum exit pressure. Hence, 
maximum isentropic compression was obtained in the supersonic 
convergent, greatly reducing the shock contribution, such that 
when combined with the pressure lift in the diffuser, the pressure 
at the ejector exit was globally higher. 


3.1.2. Constant-area mixing model 

Keenan and co-workers investigated the constant-area mixing 
ejector when they worked on the constant-pressure mixing model. 
They found that the constant-area mixing model offers better 
agreement with the experimental results than that of the constant- 
pressure mixing model. Also, the constant-area mixing model gives 
more information of the flow inside the ejector than that of the 
constant-pressure mixing model. 

In the constant-area mixing model, except for the mixing 
method and the region, the assumptions made for the analysis of 
flow within the ejector are the same as those made in the case of 
one-dimensional constant-pressure ejector flow model. 

The models proposed by Eames et al. [9] and Huang et al. [12] 
mentioned in Section 3.1.1 considered that the constant-pressure 
mixing occurred in the constant-area section. Therefore, both 
models are based on both constant-pressure mixing and constant- 
area mixing. 

Grazzini and Mariani [17] derived a theoretical model based on 
the constant-area mixing theory. In their model, the following 
hypotheses were taken: 


vapour 


Paos Tso 
Secondary vapour 


Initial interaction region 


(1) Ideal gas behavior of the superheated steam. 

(2) Same static pressure for primary and secondary flow at the 
nozzle exit plane section. 

(3) Constant-area, adiabatic mixing process. 

(4) Isentropic, one-dimensional, steady flow for both separated 
streams and for the mixed stream. 

(5) Losses are introduced by means of the isentropic efficiency for 
the primary nozzle and the diffuser. 


Using Eq. (16), the area ratio based on gas dynamic equation is: 


Ay 1 2 y-1 o 
= 14 M2 58 
Am \|\Moi Ga 2o P uan 


With a recursive method to calculate the mass, momentum and 
energy balance equations using Mach number at the end section, 
the end thermodynamic variables at the ejector exit can be 
obtained using Eq. (7). 

This model was developed for designing a jet pump device for a 
water refrigeration cycle. In order to improve the overall 
performance and the compactness of the system, they used the 
model to analyze a two-stage ejector which showed that the 
solution provides greater compression ratio while using the same 
entrainment ratio of the single stage jet pump. They also analyzed 
a jet pump with a third stage and found that the ideal gas model 
adopted could lose physical meaning. 

In order to give a more accurate prediction of the performance, 
Yapici and Ersoy [18] derived a local model based on constant-area 
mixing ejector consisting of a primary nozzle, a mixing chamber in 
cylindrical structure and a diffuser (Fig. 6). The mathematical 
description of the model is as follows: 

Primary nozzle: By using equations, such as (10), (13), (15), 
(49), and the mass flow rate at throat (Eq. (1)), they derived the 
following expressions to describe the primary nozzle pressure 
ratio (Eq. (60)) and the nozzle area ratio (Eq. (61)) 


yY 
m = pAu = PAM (59) 
RT 
P 1 1 a 
a 1 i - (60) 
po Na M (1 +¥5'M3,) 
Pee a 
y-1 = -1 72 

Api _ 1 1- mT 1+7 (61) 
Ap,nt Mp1 1--— (y-1)M, 1 + vt Me, 


2na +n (V—1)M2, 


Mixing section: For clarity, the control volume selected to 
analyze the flow in the mixing chamber is shown in Fig. 6. From 


4 


0 noe 
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vapour 


Aerodynamic throat 


Fig. 6. Ejector model applied by Yapici et al. 
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Fig. 7. Control volume selected to analyze the flow in the mixing chamber. 


isentropic gas dynamic equation, the static pressure ratio P,1/Pp is, 


Pa _ (Ta) (An) fa. Mn), ie 


where f (y,M) is the mass flow function: 


hym =m|ya +5 m] (63) 


Based on the continuity equation and energy balance for the 
selected control volume shown in Fig. 7, they derived: 


1 
Pm3 (2) GE) fa(v,Mp1) 
= + 1 64 
Pp Tpo Am3 PO, Mm) © ) (aa 
m (1407/140) (65) 
po po 
Pso — Pa f1(¥,Mp1, N) 
= 66 
Ppo Ps f3(¥,Ms1) ie 
where 
P zi vt 
Fats) = p= [1+ 45" ma] (67) 
sO 


Considering that p= P/(RT) and u? = MyRT, they derived the 
momentum equation in the x direction for the control volume 
(Eq. (68)) and the exit Mach number of the mixing chamber 


(Eq. (69)). 
Ping Am3 


Psy Ası i 2), \ 2 2 
Bo nae + yM3,) + (14 YMp1) = Pp Ap: (Lt YMims) (68) 


Mm3 = | (a? —2) + 


v-oe- 2 (69) 
where 
a= [nomo torrma (GH Ge 0 | o} 
(70) 
nonya tM (1 vow) 5 


Fig. 8 shows the control volume for the initial interaction region 
of the mixing chamber. This control volume was used for analysis 
of the flow in the supersonic regime. In order to analyze the flow in 
this region, the following additional assumptions were made: 


(1) Streams do not mix and are isentropic between sections 1-1 
and 2-2; 

(2) The secondary flow Mach number at section 2-2 is unity, 
Mg2 = 1; 

(3) The primary static pressure at the inlet is greater than that of 
the secondary, Ppi > Ps1 


Considering that Ap2 = Am3 — Ap: and Ası = Am3 — Api, together 
with Eq. (61), they derived 


1—(Ap1 /Am3) 1 
Age [Pe | 1 14%! (i y-1 T 
Ap Ap1 /Am3 Mp1 1 +% M2, Maly + 1) 
— 1)M2 y- 
' m 5 (72) 
2a + (Y — 1)M2, 


In addition, 


Thus, M,; can be found by using the above two equations. 
Considering that M,2 = 1 and Ası = Am3 — Api, they derived the 
inlet static pressure ratio: 


(Pp2/Ppo) (Ap2 /Ap* ) , 2\ 5 
Py Poi lPpo) Apt Ape) (1 yM2,) (1 + yM3; ) 


P 1-Ap Ama) [74 «afl _ PaPa Y 
p3 Ap Ana [a E YMS,) — ra ETA 


(74) 


Diffuser section: By applying the Eq. (22) to calculate the 
diffuser pressure ratio, they derived the following expression for 
the pressure ratio through the ejector 


Pao _ Pao Pm3 Ppi Poo 


75 
Ps Pm3 Ppi Poo Pso ( ) 


The authors derived the model for ejector optimization. 
Meanwhile, they carried out the parametric study by applying 
the model and derived the optimum design curves for constant- 
area ejector design. Further, they compared the analytical results 
with similar model developed by Sun and Eames [19] which was 
based on the constant-pressure mixing model, for the same 
operating temperatures of the ejector refrigeration system, and 
found the optimum coefficient of the performance and area ratio 
are greater. 


3.2. Two-phase flow model 
It should be pointed out that the above mentioned models are 


based on the assumption that the primary gas is superheated and 
the flow in the ejector is compressible single-phase flow and 
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Fig. 8. A control volume selected for analysis of flow in the initial interaction region. 


a normal shock wave may develop inside the diffuser. However, 
under many application conditions, the fluid mixture easily 
undergoes phase change and a condensation shock may develop. 
And because of the overall rise in pressure in the diffuser as well as 
the additional rise in pressure caused by the shock wave, the fluid 
mixture is easy to condensate. On the other hand, the secondary 
fluid will accelerate to sonic speed at the “hypothetical throat” 
which will reduce the pressure. Thus the sub cooled or saturated 
secondary fluid will flash into vapor and form either a two-phase 
mixture or pure vapor in the region occupied by the secondary 
fluid [20]. In addition, when the ejector works as a pump using 
high pressure liquid to entrain gas or higher pressure vapor to 
entrain liquid, under these cases, the flow is two-phase flow. 
Taking account of these cases, the models reviewed in this part are 
capable of properly simulating ejectors containing two-phase 
flow. 

Sherif et al. [20] derived an isentropic homogeneous expan- 
sion/compression model to account for phase change due to 
expansion, compression and mixing while attempting to provide 
a more accurate description of the constant-pressure mixing 
process. In this model, the primary fluid is a two-phase mixture 
and the secondary fluid is either a sub-cooled or saturated liquid 
having the same chemical composition as the primary fluid. 
Therefore, comparing with single-phase flow model, the biggest 
difference in the model is the introduction of dryness of the fluid 
in the calculation of the specific volume, enthalpy and entropy. 
For clarity, the ejector model used by Sherif et al. is shown in 
Fig. 9. 

Primary nozzle: Different from single-phase flow model, in the 
model, the primary nozzle is segmented into two parts, i.e. the 
converging section and the diverging section. For given inlet 
pressure, temperature and quality, the specific volume, enthalpy 


/ “aia 


Converging-diverging 
primary nozzle 


A Ra 


Constant pressure , 
Nozzle throat mixing chamber Ejector throat 


and entropy can be calculated by Eqs. (76)-(78). 


Uni = Unit + Xni(Vni.g — Unis) me 
hni = hni f + Xni (hni g — Ani) (27) 
Sni = Snif + Xni(Snig — Snif) ia 


Considering the continuity and energy equation in the 
converging section along with the definitions of the Mach number 
and the speed of sound, the isentropic flow requirement, and the 
equations of state, the enthalpy and density are, 


PniâniUni = PntAntUnt (79) 
1 2 1 2 
Agi += UG; = Ant + 5 Une (80) 
Snt = Sni = Snt,f + Xnt(Sntg — Snt,f) (81) 
Ant = Ante + Xot (Matg a ote) (82) 
1 1 1 1 

+ Xnt (83) 

Pot Pott Porg Pott 


By assuming the pressure at the nozzle throat Py and using 
Eq. (81) to compute the quality at the throat x,,. The quality can 
then be used to compute the specific enthalpy and density at the 
throat using Eqs. (82) and (83), respectively. And then the Eqs. (79) 
and (80) can be solved in the velocities um and unt. The speed of 
sound at the throat can be computed from Eq. (6) and be compared 
to the velocity term unt. If the two quantities are found different in 
value, another value of pressure is tried until convergence. Then 


Diffuser 


Fig. 9. Ejector model used by Sherif et al. 
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the mass flow rate of the primary fluid can be calculated by 
Eq. (59). 
Using the same method by assuming the pressure at the nozzle 
exit, the state parameters at the nozzle exit can be solved. 
Secondary flow: If the secondary flow is all liquid, Bernoulli’s 
equation can be applied between the inlet and exit of the 
secondary flow tube, 


1 
Pso = Pse + 5 Psose (84) 


Assuming that the liquid density is nearly constant throughout 
the secondary fluid path, i.e. Pse % Pso, the secondary fluid flow rate 
can be calculated using Eq. (59). 

For given entrainment ratio, Eqs. (84) and (59) can thus be 
solved to produce the pressure and velocity at the secondary flow 
path exit. 

Mixing chamber: With the assumption of constant-pressure 
mixing, the continuity, momentum, and energy equations in the 
mixing chamber are 


Mp(1 + @) = PmeAmeUme (85) 

(Phe Pse)Ane = Mp(1 t W)Ume MpUne — MpWuUse (86) 
1 2: 1 2 1 2 

(Ane j 7 Une) y (Ase t 7 Use) @ (hme j 2 Une) (1 w) (87) 


Similarly, the state parameters at the exit of the mixing 
chamber can be obtained by using the state equations similarly to 
Eqs. (81-83) as well as the equations defining the Mach number 
and the speed of sound. 

Diffuser section: The conditions at the mixing chamber exit 
were assumed identical to those at the diffuser inlet and a shock 
wave is assumed to occur, then the mass, momentum and energy 
balance in the diffuser are, 


PmeUme = PssUss (88) 
Pre + Prene = Pss + Pela (89) 
hme + i, = hss + si, (90) 
p = p(P,h) (91) 


Again, by assuming the pressure, the state parameters can be 
calculated using above equations. 

For given entrainment ratio and inlet fluid states, the model can 
be used to determine the distribution of thermodynamic states as 
well as the area of the mixing chamber exit. By analyzing the 
performance of the ejector applying the model, the authors found 
that the variation of compression ratio with entrainment ratio 
follows a monotonically decreasing trend. The more energetic 
primary fluid appears to be more suitable for inducing high 
compression ratios than high entrainment ratios. 

Cizungu et al. [21] derived a two-phase thermodynamic model 
to calculate the entrainment ratio shown in Eq. (92). 


Ahexp frst 
y Mna- (1 +32 E 
o= Sem, i ( - DEn d ) (92) 
(1+ E) -vienna 
where fm represents the friction factor, and Em is the momentum 


loss factor, the enthalpy differences are given by Ahexp = hpo — ħp1,is 
piis and Ahexs = hst — hsıis-The authors simplified the equations 


mentioned above for simple calculation: 


| Ahexp 
= —1 93 
w Ahon” (93) 


where y= Nexp/Ncom represents the quality grades of the jet 
compressor, it ranges between 0.689 and 0.81. 

This model can be used both for single-phase flowing ejector 
and two-phase flowing ejector with single-component or two- 
component as the working fluid. For given boundary conditions the 
optimum geometry of the ejector can be calculated. For given 
geometry of the ejector, the optimum thermal conditions can be 
determined. They used the model to analysis the interaction 
between the entrainment ratio and all relevant dimensional 
parameters and compared the results with experimental data. 
Good agreement was found. They also found that the dimensions of 
ejector configuration have a dominant influence in deciding the 
operating range. 

Another kind of widely used two-phase flowing ejector is that 
employs liquid as primary fluid and gas as secondary fluid or vice 
versa. This kind of ejector found widely application in industry, 
such as steam-driven jet pump used in modern thermal power 
plant to remove non-condensable gases from the condenser, 
ejectors used to entrain and pump fumes and dust-laden gases 
which otherwise are difficult to handle and so on. Therefore, 
Beithou and Aybar [22] employed a one-dimensional control 
volume method to develop a local mathematical model for the 
performance prediction of the steam-driven jet pump. In the 
model, the conservation of mass and momentum equation as well 
as the entropy equation were applied. As the fluid is uncompres- 
sible, therefore, the Bernoulli equation can be used to calculate the 
entrained water velocity at the nozzle. Regarding to the mixing 
section, the steady-state energy equation was used to calculate the 
velocity of the mixture with the assumption of adiabatic mixing 
process as well of no potential energy change. Take the head loss 
into account, the Bernoulli equation was used in the diffuser to find 
out the pressure and velocity profiles along it. The mathematical 
model is shown as follows. 

Steam nozzle: 


Conservation of mass : 


XO piuiAi = Ñ PetteAc (1) 


ðu ðu ƏP 


Conservation of momentum : P + pu K + ae 0 (94) 


The governing equations of steady-state subsonic or supersonic 
isentropic flow: 


d 

aq (PUA) = 0 (95) 
du 1 dP 

tr podi (96) 


Also the entropy equation Eq. (11) is used. 

Water nozzle: 

By using Bernoulli equation between the water tank and the 
exit of the water nozzle and the continuity equation, the water 
velocity and water flow rate can be calculated by Eqs. (97) and 
(59): 


Pr Pe ) 
ico =4 (2 97 
a (o Pwe 7) 


where, the subscript “T” defines the water properties in the supply 
water tank, and the subscript “wc” shows the water properties at 
the inlet of water nozzle. 
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Mixing section: The steady-state energy equation with 
assumption of no potential energy change is similarly to Eq. (30). 

Diffuser section: Applying Bernoulli and considering the head 
loss, they derived: 


u Py, wu P 
dod Bee Eh 98 
2 Pa 2 Pe r ( 


where h; is the head loss which is given by, 
2 
n=- G92- n] (99) 


where 7 is the pressure recovery coefficient. 

In the steam nozzle, the iterative finite difference procedure 
with under-relaxation is used to approximate the derivatives in the 
steam nozzle governing equations that are expressed in terms of 
the mesh points. 

The initial condition is the state parameters and geometrical 
parameter of the inlet section of the steam subsonic nozzle. The 
nozzle shape is such that the area variation with distance along the 
nozzle is given by a quadratic function. In the subsonic nozzle, the 
effect of the density changes is negligible. Assuming p; = pı gives 
the initialization of density at every grid point, and then calculates 
the velocity and the pressure at each grid point for the first 
iteration by, 


PiWA, 
ui = ——— 100 
' pii ( ) 
P; = Pia — pjui(ui — Ui-1) (101) 
n+1 n Pi uy n 
j= pe tr Pils) Pi (102) 


where the index 1 shows the initial condition, n shows the density 
at the previous iteration, and the under-relaxation factor r is to be 
chosen less than 1. Then the pressure, velocity and density can be 
updated by following equations: 


A 
P, =P; (2) (103) 
P; — Py 
E Ë Li 2] (104) 
i hie | 
u,A 
pr! = ph + r(e m= pt| (105) 
i adia 


The goal of the model is to evaluate the pressure, density and 
velocity distribution along the ejector. The heat transfer rate from 
steam to water was obtained using the direct-contact heat transfer 
coefficient as a function of the length of a conical steam jet plume. 
The condensation profile, non-condensed steam percentage in 
volume versus non-dimensional length of mixing section was 
obtained from the heat transfer rate profile. They compared the 
calculated pressure profile of the model with Cattadori et al.’s [23] 
experimental pressure profile which found that the experimental 
and calculated pressure distributions are in good agreement in the 
mixing nozzle and diffuser, and there are some difference in steam 
nozzle. 

The above mentioned local model models the flow in the 
injector at each point, and taking into account major phenomena 
such as condensation of droplets, heat and momentum transfers, 
viscous dissipation on the wall non-adiabatic and non-equilibrium 
flow. Such a model is more realistic, but needs a perfect knowledge 
of the previously mentioned phenomena as well as their mutual 
interactions that makes the model very complex. Take these in 


Fig. 10. Control volume for the modeling of the mixing chamber. 


consideration, Deberne et al. [24] derived a global volume model 
(Fig. 10) for steam ejector performance with the consideration of 
void fraction, area contraction ratio and condensation rate in the 
conservation of mass, momentum and energy equations based on 
steady and one-dimensional assumption. 

By defining the following parameters: 


Void fraction: & = ee (106) 
1s“1s 

(w Pithi + 1) 
Area contraction ratio: QR =- ——2—— 107 
Ais + Ay (107) 

: X1 Ly 
Condensation rate: R= | ———— (108) 

Cpt ATw P 


The conservation equations applied to the global volume lead 
to: 


1 
Mass: puz = (1 + œ) z381 Pistis (109) 
2 1\., 
Momentum: (1+@)u2 4 P+ (1 P* 
( ) á man : ( 2) ) 
Pis 
= (Us + WUjs) 4 110 
(Us + @ls) E1 Piss ( ) 
2 2 2 
Energy: (1+) hs Z| ohn aa [ms f "l (111) 
State equation: h = h(P2, p2) (12) 


where the equivalent pressure P* are obtained by the empirical 
correction numerically introduced from the experiment results at 
the CETHIL: 


P* = 1.01R?°”? x Pis (112) 


wherein, the subscript of 1s refers to mixing section inlet. The 
subscripts 11 and 2, respectively, denotes the water mixing section 
inlet (water nozzle outlet), the mixing section outlet planes. 
AT = Tsat — Tii (Tsat is saturation temperature at pressure P) is the 
inlet liquid sub cooling, Ly is the latent heat of water, C,) is the heat 
capacity of liquid and x, is the vapor quality at the mixing inlet 
section. 

Shock wave: 

They discussed the shock wave under two conditions: (i) too 
high void fraction, the fluid will not be condensed enough to reach 
a complete condensation across the shock wave and (ii) the void 
fraction is too low, the fluid velocity will be too low (lower than the 
sound speed) and no shock wave can occur. 

The first sound speed corresponds to the case where response of 
droplets is negligible (droplets are frozen with no mass and 
momentum transfers), and is similar to the sound speed in pure 
gases: 


Cp = 4/ YRgTs 


(113) 
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The other sound speed, presented in Young and Guha [25] is 
called the full equilibrium sound speed, where all processes are in 
equilibrium during shock wave: 


(114) 


where c = Cps + ((1 — €)/é)cpi, in the model, the normal shock wave 
is consider in full normal shock wave, fully dispersed with 
discontinuity in flow properties and partly dispersed shock wave. 
For shock wave, and neglect the effect of the wall on the wave, 
surface area of sections 2-2 and 3-3 are taken to be equal, and then 
the conservation equations are as follows: 


Mass: (3U3 = 2U2 With p= €p2zs + (1 -— €) py (115) 
Momentum: 33 + P3 = 0,u3 + P2 (116) 
12 
Energy: p3u3( h3 + 543 
u2 
= Po + u2[e2Pzsħhzs + (1 — £2) Paihai] (117) 
State equation: h3 = h(P3, p3) (12) 


Diffuser: By introducing a pressure loss term in the Bernoulli 
equation, then the Bernoulli equation is, 


1 2 


P44 z (3U4 


Ps 4 + p38 NaPs¥3 (118) 

Once the inlet conditions are known, as well as the steam 
injector geometry, the condensation rate R is calculated. The term 
P* is calculated with the equivalent pressure ratio t* and then 
all the other constants. If no solution for the system of Eqs. 
(116) + (12) + (117) + (118) can be found, the calculation aborts if 
the flow velocity uz is lower than sound speed upstream of the 
shock wave to ensure the assumption (two-phase supersonic 
upstream flow, shock wave and liquid downstream flow). Physical 
quantities of the flow are then calculated upstream and down- 
stream of the shock wave and the minimum void fraction condition 
is a posteriori verified (to be sure that a normal shock wave exists). 
The outlet quantities are, in this case, calculated. 


Table 1 
Features of thermodynamic model reviewed in present work. 
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The authors also set up a test rig to verify the model and 15% 
performance prediction accuracy had been found. Then a series of 
parametric investigation were carried out by applying the model to 
predict the influence of different geometrical parameters of the 
ejector (mixing section throat diameter, position of the liquid pipe, 
etc.) and the influence of physical parameters (temperature, liquid 
and steam pressure, etc.) on the ejector performance. 

The thermodynamic models summarized in this part are 
reviewed in Table 1. 


4. Dynamic model 


Despite the usefulness and the remarkable progress the 
thermodynamic models provided for the general understanding 
of ejectors, this kind of models were unable to correctly reproduce 
the flow physics locally along the ejector. It is the understanding of 
local interactions between shock waves and boundary layers, their 
influence on mixing and recompression rate, that will allow a more 
reliable and accurate design, in terms of geometry, refrigerant type 
and operation conditions. A way of achieving this objective at a 
reasonable cost is through computational fluid dynamics (CFD) 
modeling approach which provides a better understanding of the 
hydrodynamics of the ejector. Such models are called dynamic 
model (listed in Table 2). Although dynamic model needs more 
time compared with thermodynamic model, it is more related to 
the true physics of the problem. Therefore, dynamic model is likely 
to give better agreement with experiment results. With the 
development of computer science and mathematical method, 
more and more researchers have focused on this method with 
attempt to obtain more accurate information in the ejector. 


4.1. Single-phase flow model 


The single-phase flow model includes gas-gas ejector and 
liquid-liquid ejector flow model. The former considers the flow as 
compressible while the latter as incompressible. 

Because of poor mesh resolution, the early CFD models [26-28] 
are failed to track the shock wave and poor agreement with the 
experimental results. Rusly et al. [29] simulated the flow through 
an R141b ejector by using the real gas model in the commercial 
code, FLUENT. The effects of ejector geometries were investigated 


Ref. Media Gas property Multiphase Mixing section Mixing mechanism B.c.4 Les Other relations 
flow model 
Ideal gas Real gas NXP? Section y-y CPM> CAM‘ 
[9] Steam aff Compressible single-phase ,/ <i] S Tp, Pp, Mp, Ts, Ps, ii ala Eqs. (15-21) 
[10] Steam if Compressible single-phase y Uy any ee Pa w Eqs. (15), (23-28) 
[12] R141b K Compressible single-phase Si V y Tp, Pps Ts, Ps, PŁ,At Ane A1 Eqs. (6-10), (13), (29), (30) 
[13] R717 all Compressible single-phase Ni N Tı, T2, T3 Eqs. (31-41) 
[14] R134a° al Compressible single-phase N V Tog Py Uh, dP l Eqs. (42-47) 
[15] R142b al Compressible single-phase N V Inlet conditions Eqs. (48-51) 
[16] R142b aif Compressible single-phase NI al Inlet conditions [2 Eqs. (5), (6), (52-57) 
[17] Steam Compressible single-phase ./ Al ail Tay Poy Ueo Pe Eqs. (1)-(3), (7), (16), (58) 
[18] R123 al Compressible single-phase V V N Byl a Eqs. (10), (13), (15), (49), 
(59-75) 
[20] R134a all Homogeneous model al Tiny Heyl, IPs Prt, Pya Eqs. (6), (59), (76-91) 
[21] NH3 wll Two-phase I Usp Ue Eqs. (92), (93) 
[22] Steam-water Condensation Inlet conditions Eqs. (1), (11), (30), (94-105) 
[24] Steam Condensation V Inlet conditions Eqs. (12), (109-118) 


* NXP: nozzle exit plane. 

> CPM: constant-pressure mixing. 

CAM: constant-area mixing. 

B.C.: boundary conditions. 

e LC.: initial conditions. 

” Other media used are R152a, R290, R600a, R717. 


c 


d 
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Table 2 
Features of dynamic model reviewed in present work. 
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Ref. Media Gas property Multiphase flow model Turbulence model Mixing mechanism Other relations 
Ideal gas Real gas CPM CAM B.C. 

29] R141b V Compressible single-phase Realizable k-e N E Eqs. (4), (119-122) 
31] Air M Compressible single-phase RNG k-e, SST, k-w all Io Phy To Peo Pe Eqs. (4), (119-126) 
34] R142b y Compressible single-phase k-e, SST, k-w y I Po lo Pala Eqs. (4), (119-126) 
35] R141b y Compressible single-phase RNG k-e yi Ry By lP Eqs. (4), (119-124) 
36] Steam V Compressible single-phase Realizable k-€ y all PPA Js Eqs. (4), (119-124) 
37] Water-air Air ASM Standard k-e Up, Ps, Pe Eqs. (124), (127-130) 
38] Water-air Air Mixture Up, Ps, Pe Eqs. (127-130) 

39] Water-air Air ASM Eulerian Standard k-e URBE Eqs. (127-133) 

41] Air-water Air mixture Standard k-g Mp, Ms, Pe Eqs. (127-130) 

42] Water-vapor Separated two-phase Pp, Tp, Ps, X Eqs. (134-149) 


numerically. The CFD’s results were validated with experimental 
data provided by others [12] and good agreements were found. 

As it is often the case in transonic compressible flows involving 
shock reflections and shock-mixing layer interaction, the choice of 
the turbulence model and grid refinements are critical points [30]. 
Therefore, Bartosiewicz et al. [31] used the CFD model aiming at 
validating the choice of a turbulence model for the computation of 
supersonic ejectors in refrigeration applications. In order to reduce 
the complexity of the model and to be able to use the available 
experimental data, they chose air as the working fluid. Six 
turbulence models, named k-e, realizable k-e, renormalization- 
group (RNG) k-e, RSM and shear-stress transport (SST) k-w have 
been tested and compared with measurements of Desevaux and 
Aeschbacher [32] and Desevaux et al. [33]. The compressible 
steady-state axisymmetric form of the fluid flow governing 
equations are shown as Eqs. (119-122) together with ideal gas 
law Eq. (123). 


A204? (puiA) =0 


Conservation of mass : JE t Ox, (119) 
Conservation of momentum : E (pui) + i (puiu) 
: Jt pui Ox; PUjU; 
oP ƏTij 
= -at Ox, (120) 
. ð ð 
Conservation of energy : ap (PE) 5 (u (pE + P)) 
Xi 
3 oT. = 
= V (Aeff a) + V(Uj(Tij)) (121) 
OX; 
with 
E Ou; du; 2 OU, 
Tij nor(k H 7) 3 Mett OX, bij (122) 
P 
P= (123) 


As the k-g, RNG k-e and k-w models are based on the Boussinesq 
hypothesis which makes these approaches need relatively low 
computational cost associated with the determination of the 
turbulent viscosity. However, the main drawback of the approach 
is the assumption that the turbulence is isotropic. While the 
Reynolds stress model (RSM) does not rely on the Boussinesq 
assumption, but the associated CPU cost is relatively high. 
Comparing the vast simulations by applying these approaches, 
Bartosiewicz and his coworkers found that the RNG k-e and SST 


k-w models were very promising for ejector analysis. Therefore, 
they compared these two approaches in detail. 
The RNG-k-€ model (RNG) 


a a ðe E e2 
OX} (peui) OX; ete a T Cie Kok Cas P k Re (124) 
With R, =H" C—ata) y= Sk/e, no=438, B= 0.012, 
Ci, = 1.42, G, = 1.68, CG, = 0.0845, and w= pC,.(k?/e). 
The standard and SST k-w models 
ð a) POE _ Mr) 90] 
ge (0) + (pou) =g (H+ BE) E| + Gy Yu (128) 
With u, = & bk and 
ð _ 9 o 9 _ He) 00] | 
gp (0) + g (POM) = ae (w+) al (oat 
-= F1)Do (126) 


The investigation has shown that the RNG k-e and SST k-w 
models were the best suited to predict the shock phase, strength, 
and the mean line of pressure recovery with the SST k-w model 
shown better performances in term of stream mixing. 

Then they used the validated model to reproduce the different 
operation modes of a supersonic ejector, ranging from on-design 
point to off-design to capture the shock-boundary layer interac- 
tions, such as boundary layer separation, flow separation and 
recirculation at the motive nozzle exit and the diffuser. The results 
turned out that CFD is an efficient diagnosis tool of ejector analysis 
(mixing, flow separation), for design, and performance optimiza- 
tion (optimum entrainment and recompression ratios). 

In Bartosiewicz et al.’s [34] work, they improved the previously 
mentioned model by using real gas thermodynamic and trans- 
portation properties derived from NIST-REFPROP database to 
perform simulation studies on ejector operation with R142b. They 
observed over-expanded and a strong oblique shock at the 
primary nozzle exit which induces a mixing layer separation 
associated with a corresponding energy loss. Consequently, this 
separation creates an adverse pressure gradient at the secondary 
inlet, giving rise to favorable conditions for an outflow. The one- 
dimensional model cannot take into account such irreversible 
losses induced by oblique shock wave reflection and flow 
separation. This local CFD modeling takes into account shock- 
boundary layer interactions in a real refrigerant. The numerical 
results contribute to understanding the local structure of the flow 
and demonstrate the critical role of the secondary nozzle for the 
mixing rate performance. 
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Zhu et al. [35] employed the CFD technique with RNG k-e model 
to investigate the effects of two important ejector geometry 
parameters: the primary Nozzle Exit Position (NXP) and the mixing 
section converging angle 0, on its performance. The flow inside the 
ejector is governed by the compressible steady-state turbulent 
form of the flow governing equations. They found that the 
optimum NXP is not only proportional to the mixing section throat 
diameter, but also increases as the primary flow pressure rises. On 
the other hand, the ejector performance is very sensitive to 0 
especially near the optimum working point. The entrainment ratio 
varies as much as 26.6% by changing 8. In order to better maximize 
the ejector performance, a relatively bigger@ is required when the 
primary flow pressure rises. These findings are good for guiding the 
adjustment of NXP and 8 in order to obtain the best ejector system 
performance when the operating conditions are different from the 
on-design conditions. 

The previously mentioned dynamic models are based on 
constant-pressure mixing ejector. Pianthong [36] employed the 
CFD with realizable k-e turbulence model to predict the flow 
phenomena and performance of both constant-pressure mixing 
and constant-area mixing steam ejectors with application in 
refrigeration system. The results indicate that CFD can predict 
ejector performance very well and reveal the effect of operating 
conditions on the effective area that is directly related its 
performance. Besides, they found that the flow pattern does not 
depend much on the suction zone because the results of 
axisymmetric and 3D simulation are similar. 


4.2. Two-phase flow model 


There are two approaches for the numerical calculation of 
multiphase flows: the Euler-Lagrange approach and the Euler- 
Euler approach. In the latter approach, different phases are treated 
as interpenetrating continua and this approach has been widely 
adopted in ejector modeling. The conservation equations have 
similar form for all phases in this approach. There are two Euler- 
Euler multiphase models for ejector modeling: the mixture model 
and the Eulerian model. 


4.2.1. Mixture model 

This model solves the momentum, continuity, and energy 
equations for the mixture while the volume fraction equations are 
solved only for the dispersed phases. The model uses algebraic 
expressions for the relative velocities and the equations for 
multiphase flows are derived by mass-weighted averaging or 
Favre-averaging method. 

The continuity and the momentum equations for mixture flow 
are, 


a 
ge Cmm) + V(QmPmuUm) = I’m (127) 
a 
Je (PmUm) + V (PmUmUm) = —VP+V(tm+ 1Im) + VtDm 
+ Pm&+Mm (128) 


The diffusion stress term VtDm represents the momentum 
diffusion due to the relative motions of the two phases while Mm is 
the influence of the surface tension force on the mixture and 
depends on the geometry of the interface. 

In order to calculate the relative velocity, it is necessary to 
calculate the average interfacial momentum source for the 
dispersed phase Mp. By neglecting the surface tension forces 
(Mm = 0) and assuming that the phase pressures are equal and 


rearranging, the equation for dispersed phase M, is 


Mp = fanto pm) E + (UmV)Um -s|} 


Um, 


a 
4 {avp oe (uy, V)um,| } 
+ {ppp [(UmV)Um, + (Um, V)Um |} 
+ {-V [@p(Tp + Tr,)] } + @pV(Tm + Tm + Tm) 


(129) 


The second term from Eq. (128) is omitted by assuming local 
equilibrium (the bubbles are rapidly accelerated to the terminal 
rise velocity). The third term corresponds, in rotational motion, to 
the Coriolis force while the last term is the contributions of the 
viscous, turbulent and diffusion stresses. These are usually 
neglected. 

By neglecting the virtual mass effect, lift-force effect, basset (or 
history) effect except the viscous drag and assuming spherical 
bubbles, the momentum transfer induced by the fluid-fluid 
interaction force is 


1 
Mp = — 3 PApCp|Ucp|Uep 


g + (UmV)Um + = (130) 


= Vo(Pp Pm) ot 


a 
where Uep is the relative velocity, Ucp = Uc — Up. 

Kandakure et al. [37] used the algebraic slip mixture (ASM) 
model to study the hydrodynamic characteristics of the ejector. 
The model accounts for the relative velocity between the 
continuous and dispersed phases as the interfacial momentum 
exchange term. In the simulation, water was taken as the primary 
fluid and air as the entrained fluid. The air phase was assumed to 
obey the ideal gas law. Standard k-e model was used for modeling 
the turbulent behavior of the flow. The slip between the phases 
was taken as some percentage of local primary velocity. They found 
that at low value of area ratio (ratio of the mixing chamber’s throat 
area to nozzle throat area), due to the larger diameter of the water 
jet, the annular area available for air flow reduces, causing 
recirculation of the entrained air within the converging section of 
the ejector. On the other hand, for higher values of area ratio, due to 
smaller diameter of the water jet, the momentum transfer to the 
air decreases and all the entrained air cannot be forced through the 
throat. As a result, the net air flow rate going into the throat for 
both area ratios is small. Thus there is an optimum area ratio for 
the maximum air entrainment rate. The air entrainment rate 
correlates with pressure difference between the air entry and 
throat exit for a wide variety of ejector geometries and operating 
conditions. The overall head loss factor and the ejector efficiency 
can be predicted a priori. 

Kim et al. [38] also investigated the hydrodynamic character- 
istics of the liquid-gas (water-air) ejectors with mixture model 
and experiments focusing on the gas hold-up performance. They 
found the gas phase hold up increases with increasing liquid 
circulating rate but it decreases with increasing liquid level in the 
column. The gas suction rate increases with increasing the liquid 
circulating rate but it decreases with increasing the liquid level in 
the column and nozzle diameter. 

In order to get maximum entrainment of the secondary fluid, 
Yadav and Patwardhan [39] used the ASM model to determine the 
optimum geometry of suction chamber of the ejector which is 
worked as a contactor for two fluids mixing. It is found 
experimentally that the entrainment of the secondary fluid was 
more when straight portion of the throat is of zero length [40]. 
Hence, the geometrical model of the ejector adopted by Yadav and 
Patwardhan is shown in Fig. 11 with zero length of throat. Also, 


1776 S. He et al./Renewable and Sustainable Energy Reviews 13 (2009) 1760-1780 


(a) Inlet for primary fluid (water) 
Nozzle 


Inlet for secondary fluid (air) 


Suction chamber 


‘Axis of symmetry 


Diverging section 


Doutiet 


Outlet 


(b) Velocity inlet 
Pressure inlet 


Location for calculating 
area available for flow of 
secondary fluid 


Projection ration=L7y/Dr 
Area Ratio=((Ds’-Dyx" Dy) 
Angle of converging section=0 


Pressure outlet 


Fig. 11. Geometrical model adopted by Yadav and Patwardhan. 


water was taken as the primary fluid and air as the secondary fluid 
with air assumed to obey the ideal gas law. The standard k-e 
turbulence model (per phase) was used to calculate the turbulence. 
Then they carried out a large number of CFD simulations using the 
model to understand the effects of the geometrical parameters of 
the suction chamber on entrainment of the secondary fluid. They 
found that the geometry of the suction chamber has significant 
effect on the rate of entrainment of the secondary fluid. After 
studying the three parameters: the projection ratio (PR: distance 
between nozzle tip and entry to throat), area ratio and angle of the 
converging section based on the model, Yadav and Patwardhan 
recommended that a PR of five produces the highest rate of 
entrainment of the secondary fluid. Increasing the PR beyond five 
does not lead to any significant improvement in the suction 
capacity of the jet. In addition, they found a suction chamber with 
(D2 — D2) /D2 = 6.6 (D,: diameter of suction chamber, Dp: diameter 
of nozzle exit) gives the largest entrainment. A larger suction 
chamber leads to a reduction in suction capacity. The angle of 
converging section is recommended in the range of 5°-15°. 

The previously mentioned models are taken water as primary 
fluid and air as secondary fluid. Balamurugan et al. [41] used the 


mixture model with standard k- turbulence model to investigate 
the hydrodynamic characteristics of ejectors using air as the 
primary fluid and water as the entrained fluid. Both air and water 
were considered at room temperature. The air was assumed to 
obey an ideal gas law. They observed that the liquid entrainment 
rate increases with an increase in the liquid level and the gas 
velocity. The entrainment rate was found to be the highest 
corresponding to area ratio of about 0.155. 


4.2.2. Eulerian model 

The Eulerian model solves the momentum and continuity 
equations for each phase separately. The coupling of the equations 
is achieved through the pressure and inter-phase exchange 
coefficients. Momentum exchange between the phases is also 
dependent upon the type of mixture under consideration. 

The continuity equation and the momentum equations for the 
Eulerian model are the same as the mixture model. The average 
interfacial momentum source for the phase k is 


n 
Mk = So(Ri + IM jRU je) + Oe Px (FR + Fiif + Fum) (131) 
Al 


where Fẹ is an external body force, Fiifex is a lift force, Fym,x is a 
virtual mass force, Rj, is an interaction force between phase i and 
phase k and p is the pressure shared by all the phases. The lift force 
acts on a dispersed phase due to velocity gradient in the 
continuous-phase flow field, is: 
Fiiftp = —0.5 0.Op(Uc — Up) x (V x Uc) (132) 

For multiphase flows, the virtual mass effect occurs when a 
dispersed phase accelerates relative to the continuous phase. The 
inertia of the continuous phase mass encountered while accel- 
erating dispersed phase exerts a “virtual mass force” on the 
dispersed phase, which is given by: 


D 
Fym,p = 0.50 Pe(RE (Uc = Up)) (133) 


Eulerian model assumes dispersion type flow where the 
secondary fluid is distributed in the primary fluid. However, in 
Yadav and Patwardhan’s [39] work, the ejector was arranged such 
that the water jet was directed downwards. The entrained air 
flowed around the water jet in the annular space between the 
water jet and the ejector. They found no bubble formation inside 
the ejector and both the fluids flowed co-currently. Since there was 
no bubble, the virtual mass, which arises from acceleration, and lift 
force, which arises from unbalanced pressure forces around the 
bubble, do not exist. Hence, the Eulerian model is not a suitable 
choice for studying the coaxial type of flows in ejector. 

The previously mentioned models consider the primary and 
secondary fluids are of different chemical compositions with non- 
condensable gas. However, when both fluids are of the same 
chemical composition, the vapor is condensable in the fluid, there 
may be need a different model to simulate the hydrodynamic 
performance of the ejector. Therefore, Narabayashi et al. [42] 
derived a separate two-phase flow model to analysis the 
performance of a steam injector entrained by high pressure water 
in order to check the feasibility of a large-scale steam injector for 
which a demonstration test was not able to be conducted. In the 
model, attention was paid on the radial heat transport enhance- 
ment by the turbulence because of relatively large-scale water jet 
being used. Also, the interface activation was taken into account. 
The model including steam phase, water phase and the interface is 
shown as follows. 


S. He et al./Renewable and Sustainable Energy Reviews 13 (2009) 1760-1780 1777 


Steam: 
div(p,Ug) = -mg (134) 
: > oP 
div(p,UgVg — uggradVg) = -F (135) 
; = oP 
div(pzUgWg — uggradWg) = -z" Ci|W, — Wi| (Wg 
— Wi) 136 
; > DP 
div(p,Ughg — I'hggradhg) = -DTT p 137 
div(p,Ugk, — kggradkg) = Gkg — pgég 138 
F = E 
div(p,Ugég — l'eggrad eg) = Cy A — C2 PpE2/kg 139 
Water 
div(@ Ü) = —mg (140) 
: = oP 
div(p,U,;V, — PujgradV,) = =ar (141) 
: = oP 
div(p,U,W, — CujgradW;) = =a. C\|We — Wi| (We 
— Wi) + mg We 142 
div(e,U,h; — Mhygradh,) = 0 143 
div(p,0,k, — Ikgradk;) = Gk, — p,& 144 
div(p,U,e, — Pegrad £1) = C4 jë Cki — Copier /ki 145 
1 


4.2.3. Interface 

Assuming direct condensation onto the surface of the water jet, 
the heat transfer coefficients of direct condensation are given for 
each interface cell surface, by using the ideal condensation based 
on the Clausius-Clapeyron equation and Bernoulli’s theorem: 


Ha = VOR? /((Ug(Ug — U1))(Tsat(Tsat — Ti))) (146) 
Mg, = Ha(Ts — Tyi)Ai/hig (147) 
div(o Üh, — Chigradh,) = mghg (148) 


The interface velocities between the phases for each interface 
cell surface are: 


Wii = (TgWg + TW))/(Ty + T1) (149) 


where Tg is the gas shear stress and T, is the liquid shear stress, C; is 
the shear coefficient of interface between gas and liquid phases 
and the subscript and i is interface between gas and liquid phases. 

They verified the model by using high pressure steam test data 
(high pressure small size model), as well as the visualized data 
such as temperature and pressure distribution data. Good 
agreements were found between the test results and the analysis 
results. Then they used the model to analysis the characteristics of 
the large-scale steam injector which shown that the steam injector 
driven jet pump could work in the high pressure range over 7 MPa, 
and discharged over 12 MPa even at the large-scale rated flow rate 
of 220 ton per hour. 


5. Empirical/semi-empirical model 


The physical phenomena involve supersonic flow, shock 
interactions, and turbulent mixing of two streams inside the 
ejector enclosure which are difficult to be modeled precisely. Also, 
the solving process of the models based on the physical 
phenomena is complicated and time-consuming. Hence, in order 
to the evaluate the ejector performance and to instruct design, 
Huang et al. [43] derived empirical correlations for ejector design 
using R141b as the working fluid after experimented on 15 sets of 
ejectors: 


hs (jm) +14300(2) 
= — _09.0517(—@ 1.4362(—@) — 4.1734 150 
Ans Ang) Re oom 
Am 2 2 2 
Z, 7P + bire + bate? + barg + b4rcrg + bsr? Tg + berg 
nt 
+ byrerg? + bgre?rg? (151) 


where r: = Pt /Pe;rg = P;/Pe and the bo ~ bg are coefficients. 

Cizungu et al. [21] worked out following correlation based on 
the experimental investigations of Dorantés and Lallemand [44] to 
estimate the entrainment ratio of the ejector: 


2.12 
ie 332(3 ap) 


Yy ép 
where 2.5<&<4, Y <6. 
However, the empirical model cannot be extent to other fluid 
and ejector. In order to derive a relationship working for more 
fluids and physical configurations, Chou et al. [45] regressed a 
relationship for the maximum flow ratio against the geometry 
based on the experimental results from literatures [9,12,46-48]: 


(152) 


Ms (E) 1 
= KyK;Ks@ (153) 
Mp Pg (Tata p s 
where K, accounts for the primary nozzle performance, 
Pe\* Pe 
Kp (F) + AFPs (z) (154) 


wherein A is the correction constant for the primary nozzle 
working range and FP; is the failure pressure ratio of the primary 
nozzle. Kç is the physical effective area ratio proposed by Solokov 
and Zinger [49], Ke =u ff - where f f = P is the ejector 
geometric design area ratio for the ejector, App is the area of 
primary flow at choking location and m is a correction constant of 
between 1.35 and 1.5; Ks accounts for the “blockage effect” of the 
ejector configuration, affecting the entrainment and mixing 
process and ¢ is an application specific correction factor which 
can be applied to an ejector operating with a given refrigerant and 
system configuration. Ks = iar, where L,/A is the flow 
aspect ratio (L, is the length of the suction chamber, i.e. convergent 
section of mixture) and A is the radius of A3); b is the Mach number 
limitation constant, 4 for R113 and R141b and 5.2 for steam, and M 
is the design exit Mach number of the primary nozzle. 

More investigators have performed detailed experiments and 
developed numerous empirical correlations to predict the entrain- 
ment ratio, gas hold-up, mass transfer coefficient and interfacial 
area based on the flow direction: vertical up-flow, vertical down- 
flow and horizontal flow. Balamurugan et al. [50] summarized in 
detail of these correlations which will not be covered here. As most 
of the models are based on liquid as the primary fluid and gas as the 
entrained fluid, Balamurugan et al. [50] developed a semi- 
empirical model to predict the liquid entrainment rate (gaseous 
as primary fluid) taking into account: the compressible nature of 
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air in the nozzle, pressure drop for two-phase flow and the losses 
due to change in cross sectional area. 


The entrainment ratio: w 
2 
= 751 (2) -9.93 (32) 2.31 (155) 


The previously mentioned empirical or semi-empirical models 
are regressed from measurements. They are limited to the range 
over which it was developed, which limits their use in investigat- 
ing the performance of new ejector fluids, designs or operating 
conditions. Hence, Dessouky et al. [5] developed a semi-empirical 
model based on the basic model similar as [10] and a large 
database extracted from several ejector manufacturers and a 
number of experimental literature studies with attempting to 
eliminate the need for iterative procedures and simple to use, also 
give more flexibility in ejector design and performance evaluation 
comparing with empirical model. 

The correlations for steam ejector performance evaluation are: 


w = aErop,¢p,4 E+ fPp*) 


(h + iP.) Se 


P (e + fln(P D) 
o = aE PeP o RPA 


wherein, a-i are the coefficients. The former is used for the 
entrainment ratio of choked flow or compression ratios above 1.8 
while the latter is used for the entrainment ratio of un-choked flow 
with compression ratios below 1.8. They also worked out 
correlations used for design: 

For choked flow: 


(157) 


Pnp = 0.13P9-73 p973 E 
p E E 
Ant/Am = 0.34P1-09 p7 1-12 970-16 Hn 
Anp/Ant = 1104p, 982 P986 -0-12 nibi 
For un-choked flow: 
Pnp = 1.02P; 0-000762 p0.99 ae 
Ant/Am = 0.32PL' Par we 
Anp /Ant = 1.22P70-81 p98! 70.0739 a 
P 7 o 


They compared the results of the semi-empirical model with 
the experimental results derived from other literatures [6,9,48,51- 
55] and good agreements were found. 

Besides, they derived an empirical model for ejector perfor- 
mance analysis [56] based on the data and method presented by 
Power [11]: 


sowa P,\°215 /PCF 
f (P;)1% \Ps TCF 


(164) 


PCF is the primary steam pressure correction factor and TCF is 
the entrained vapor temperature correction factor. They are 
defined as: 


PCF = 3 x 1077 (Pp)? — 0.0009(P,) + 1.6101 (165) 


TCF = 2 x 10 — 8(T;)? — 0.0006T, + 1.0047 (166) 


The model equations are valid only for an ejector operating with 
steam as the primary fluid and the entrained fluid is water vapor. 


These equations are valid over the following ranges: 4500 > 
Ts > 10°C, 3500 > Pp > 100 kPa, and Pexit/Ps > 1.81. 


6. Conclusion 


Various mathematical models have been proposed to assess the 
performance of ejector for different operation and design 
conditions. Generally, the mathematical models are based on 
the flow and mixing phenomena inside the system. They can be 
divided into two main categories in terms of thermodynamic 
model and dynamic model. Within the two categories, the models 
can be subdivided based on flowing phases. Both thermodynamic 
model and dynamic model have two types: single-phase flow and 
two-phase flow. The single-phase flow thermodynamic model can 
be classified into two sub-categories according to the mixing 
mechanism: constant-pressure mixing model and constant-area 
mixing model. The two-phase dynamic model can be further 
subdivided based on the calculation methods: mixture model and 
Eulerian model. Also, using the measurements data, several 
empirical/semi-empirical models to evaluate the performance of 
ejector and to guide ejector design are reviewed in this paper. 

Thermodynamic models are usually expressed in explicit 
algebraic equations and are based on the steady-state one- 
dimensional model, focused on the pressure change caused by 
supersonic shock. In these models, the detailed local interactions 
between shock waves and boundary layers, their influence on 
mixing and recompression rate are not taken into account. Lumped 
method is applied in the governing equations which lead to simple 
model with lower accuracy. In order to offer good accuracy in 
predicting the performance of the ejector by the thermodynamic 
model, some isentropic coefficients accounting for the friction loss 
were applied in the model and were determined by experiment. In 
these models, certain knowledge on the choking, shock and mixing 
should be obtained in advance and assumptions should be made 
for model simplification. The more practical models known as 
constant-pressure mixing model and constant-area mixing model 
made the model feasible. The constant-pressure mixing ejector 
provides better performance than the constant-area mixing 
ejector. While the constant-area mixing model gives more accurate 
performance prediction. In order to give more accurate prediction, 
the two-phase flowing model is developed considering the 
condensation of the secondary flow or the mixture of two streams 
in different phase states. In these models, quality is introduced in 
the calculation of the state parameters. Though the governing 
equations are basically the same as single-phase flowing model, 
the accuracy is improved. 

The dynamic model, on the contrary, accounts the turbulence 
interaction between the primary and secondary stream, the shock 
reflection and chocking. It is more related to the actual process 
occurred in the ejector and the effect of the geometrical 
parameters and operation parameters can be well explained. 
The precision of this model is thus greatly improved. 

The flow in the ejector is from supersonic to sonic, and then to 
subsonic. Therefore, the choice of the turbulence model is 
important. For compressible flowing model, it was found RNG k- 
g and k-w-SST models were the best suited to predict the shock 
phase, strength, and the mean line of pressure recovery [31]. For 
incompressible flow, the standard k- model and realizable k-e 
model were widely used for time saving. For two-phase flowing 
model, the mixture model can give reasonable results. 

It can be seen that the pressure inlet and pressure outlet 
boundary conditions are widely used in all the models whenever 
possible. Other boundary conditions such as mass flow rate or 
velocity are seldom used in models. For thermodynamic model, the 
auxiliary relations such as gas dynamic equations, the defining of 
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Mach number and sonic velocity, the state expression have to be 
employed to complete the mathematical model. 

Due to the complex nature of the partial differential equations, 
it is necessary to solve the mathematical model using numerical 
methods. Finite difference method is recognized as the most 
accurate and most universal solution technique and is widely 
employed in ejector modeling. 

Besides, model validation is an important step in model 
development since it offers the possibility of comparing simulation 
results with actual system behavior. Experiments are mostly used 
to validate the mathematical model. Besides, comparing with the 
previous numerical results is also a good validation method. 

Though a large amount of works have been conducted on 
modeling and analyzing ejector, further efforts are still needed: 


(1) To study the influence of variable isentropic coefficients which 
are taken as constant in almost all existing thermodynamic 
models. 

(2) To improve the accuracy of the model based on turbulence 
model, since by now, there still no uniform model being used. 

(3) To construct a simulation package of the whole ejector-based 
system by combining the model of the ejector and other 
components in the system. 
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